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I. INTRODUCTION 

In this article, we will be interested in non-integrability proofs of meromorphic homoge- 
neous potentials of degree —1 in the plane, and in particular in nongeneric cases. Writing 
our potential V in polar coordinates, and making the Fourier expansion in the angle gives 
us 

oo 

V(r,6)=r- 1 £ a k e lk0 . (1) 

k=— oo 

This type of potential covers many physical problems in celestial mechanics and n-body 
problems, in particular the anisotropic Kepler problem, the isosceles 3-body problem, the 
colinear 3-body problem, the symmetric 4-body problem and so on. Moreover, for such a 
potential there are strong integrability conditions, thanks to the Morales-Ramis theory 1 and 
to a very effective criterion of Yoshida 2 . Still, for such a general potential, this criterion 
will not be sufficient. This is not particularly because this class of potentials is large, but 
because there are nongeneric, very resistant cases inside. For example, if we want to study 
the integrability of V(r, 9) = r -1 /i(exp(z#)) with a polynomial h, we have a priori a potential 
with deg h + 1 complex parameters, and Yoshida's integrability criterion will restrict this 
family to a family with deg h — 1 integer parameters. Still one would like to have a finite 
list of possible integrable potentials, so as to be able to check the existence of first integrals 
one by one. Here we will present a stronger criterion in Theorems 2 and 3 which is able to 
deal with such families, and which therefore is capable to settle any integrability question 
on finite dimensional families of type (1). As an application of our method, we will apply 
this criterion in the case V(r, 6) = r~ 1 h(exp(i6)) with h e C[z], deg ft. < 3. To do precise 
statements, let us now begin with some definitions concerning homogeneous potentials and 
integrability. 

Definition 1. We consider the algebraic variety S = {(qi,q 2 ,r) £ C 3 , r 2 = q\ + g|} and 
the derivations for a function f on S 

Qf df 

— = d 1 f + r^qtdaf, — = d 2 f + r~ l q 2 d 3 f 

oqi dq 2 

where d\ is the derivative according to the i-th variable (the variables of f are qi,q2,r in this 
order). This defines a symplectic form on C 2 x S on which we consider a Hamiltonian H 
of the form 

H(pi,P2,qi,qn,r) = ^(pj + pi) ~ V(q 1 ,q 2 ,r) 
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with the associated system of differential equations 

i d d 

r = r (gig! + q 2 q 2 ), ft = j— H, pi = - i = l,2. (2) 

The potential V is assumed to be meromorphic on S and to have the following form in polar 
coordinates: 

V(r, 6) = -U(6), rcos$ — qi, rsin6 = q 2 
r 

This implies that V is homogeneous of degree — 1. We say that I is a meromorphic first 
integral of H , if I is a meromorphic function on C 2 x S such that 

{ H,J } J^[dpi H dq/ dq { H dpi 1 ) °' 

Obviously, the Hamiltonian H itself is a first integral. We will say that V is meromorphically 
integrable if it possesses an additional meromorphic first integral which is independent almost 
everywhere from H. 

Definition 2. We call c = (ci, c 2 , C3) G S a Darboux point of V if 

d d 

— — V(c) = aci and — — V(c) = ac 2 (3) 
dqi dq 2 

where a G C is called the multiplicator. Because V has singularity at C3 = ; we will always 
assume that C3 7^ 0. Because of homogeneity, we can always choose a — or a — — 1. We 
say that c is non- degenerate fa^O. To the Darboux point c we associate a homothetic 
orbit given by 

r(t) = c 3 0(t), ft(t) = Q0(t), p l {t) = c4it) (^ = 1,2), (4) 
with (ft satisfying the following differential equation 

\^-~W) +E ' E€C - 

In the following, we will often omit the last component of a Darboux point c G S as it is 
defined up to a sign (and the choice of sign does not matter) by the two first components. 

Definition 3. The first order variational equation of H near a homothetic orbit is given by 

x{t) = -LvV(c)x(f) 
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where V 2 V(c) is the Hessian ofV (according to derivations in q). After diagonalization (if 
possible) and the change of variable <p(t) — > t, the equation simplifies to 

2t 2 (Et + l)Xi - tXi = \X U 

where the A, are the eigenvalues of the Hessian of V evaluated at the Darboux point c, i.e., 
Aj e Sp (W(c)). 

Theorem 1. (Morales, Ramis, Yoshida 231A ) If V is meromorphically integrable, then the 
neutral component of the Galois group of the variational equation near a homothetic orbit 
with E ^ is abelian at all orders. If we fix the multiplicator of the associated Darboux 
point to —1, the Galois group of the first order variational equation has an abelian neutral 
component if and only if 

Sp (W(c)) C [\{k - l)(k + 2) : ken). 

If the multiplicator of the Darboux point is 0, the Galois group of the first order variational 
equation has an abelian neutral component if and only if 

Sp ( V V(c)) C {0} . 

In fact, this is not exactly the same statement as the original theorem because we allow 
r to appear in the potential and in the first integrals. 

Proof. Let Y C C 2 x S denote the curve defined by equation (4) without the singular point 
(<?i> Q.2-, r ) — (0, 0, 0), and M an open neighbourhood of V in C 2 x5 such that H is holomorphic 
on M. The Hamiltonian H is then well defined and holomorphic on a symplectic manifold 
M and the additional first integral is meromorphic on M. Hence, using the main theorem 
of 3 , the neutral component of the Galois group of the variational equation near V is abelian 
at all orders over the base field of meromorphic functions on V. The variational equation is a 
hypergeometric equation. In 6 , Kimura classifies Galois groups of hypergeometric equations 
over the base field C(t). We can use this classification as the Galois group over the base field 
C(t) is the same as over the base field of meromorphic functions because the hypergeometric 
equation is a Fuchsian equation (see page 73 of 1 ). This produces the condition on the 
spectrum of V 2 V(c). The case of a degenerate Darboux point leads to the variational 
equation 

X = xr 3 x 
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which is a Bessel equation (after a change of variables). Its Galois group over the field of 
meromorphic functions in t has not an abelian identity component except if A = 0. □ 

Note that in the case of a degenerate Darboux point, we explicitly need that the first 
integral is meromorphic including r = 0, as the variational equation is not regular singular 
at this point. The integrability condition for a non-degenerate Darboux point also holds 
for a potential meromorphic only on S* = S \ {r = 0} and meromorphic first integrals on 
C 2 x <S*. 

II. MAIN RESULTS 

In this section, we are going to state the main theorems of this article. The remaining 
parts of this paper are dedicated to their proofs. 

Theorem 2. Let V be a homogeneous potential of degree —1 in the plane. We suppose 
that c = (1, 0) is a Darboux point of V with multiplicator — 1. If the variational equation is 
integrable at order 3, then the following conditions are fulfilled 

Sp (W(c)) = {2, |(p - 1)0 + 2)} for some peN. 

If p is even then 

( d 3 V \ 2 (d 3 V\ 2 fd 4 V\ 

and if p is odd then 

d 3 V ( d 3 V \ 2 fd 4 V\ 

where the functions /1, /2, fs satisfy explicit P- finite recurrences, i.e., linear recurrences with 
polynomial coefficients. 

This theorem is a generalization of the criterion given by Yoshida for homogeneous po- 
tentials in the case of degree —1 and dimension 2. A similar theorem could be proven in 
higher dimensions, but the main problem is that Theorem 2 is almost inapplicable in this 
form. In most cases, it is necessary to study more closely the expression of the functions 
h{p)i fzip), h{p) to apply it, and for the moment, because of limitations of computing 
power, it seems only possible to do in dimension 2 (for which the computations are already 
tedious). 
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Theorem 3. The functions f\(2n), /2(2n), /s(2n) can be written as 

1511011 1511011 31731231 \ 

.67108864^ ~ 134217728™ 3 + 4294967296n 4 / 



/ 1 (2n) = e 1 (n) 



/ 22665165 22665165 298125 \ 

M ™J - e 2 (rij ^ 10 7374 18 24n 4 ~ 1073741824n 5 + 4194304nV 

/ 1740684681 1740684681 2400813907 \ 

h{ n) - e 3 {n) \ - 68719476736n 2 + I37438953472n 3 ~~ 68719476736nV 



h(n) - 1| < 10~ 5 > 100. 
With this, we can apply Theorem 2 to some concrete examples: 
Theorem 4. Let V be a potential in the plane expressed in polar coordinates by 

V(r, 9) = r" 1 (a + be ie + ce 2id + de m ) . (6) 

If V is meromorphically integrable, then V belongs to one of the following families 

V = r~ 1 a, V = r- l (a + be ie ) ) V = r~ l (ae ie + be m ) 

V = r" 1 (a + be 2ld ) , V = r" 1 (a + be m ) , V = r~ l (a + be 

with a, b G C. 



(7) 



The first three families have already known additional first integrals 7 , polynomial of 
degree 1 or 2 in p. The status of the last three families is unknown. This is not due to 
an incomplete application of the Morales-Ramis Theorem, but linked to the fact that either 
they do not possess any Darboux points, or in the last case the only Darboux point is very 
degenerate and therefore the Morales-Ramis Theorem gives no integrability constraints at 
any order, as proven in 5 . 

In practical problems like Theorem 4, studying integrability only using the Morales-Ramis 
criterion is impossible because of two facts. First we need a Darboux point of our problem; 
if we do not have any, the only thing we can do is to try to find an additional first integral 
using the direct method of Hietarinta 7 . 

The second problem is the following scenario: inside the family of potentials given by 
Theorem 4, there exist submanifolds in the space of parameters for which the potential 
possesses only one Darboux point and the eigenvalue at this Darboux point can be arbitrarily 
high. In this case, the higher variational method is required. But the constraint at order 2 
does not give sufficient conditions to conclude, and it is necessary to go to order 3. 
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But the expression of this constraint cannot be written explicitly for all possible eigenval- 
ues, only for a finite number of them. To apply this third-order criterion, we derive P-finite 
recurrences and asymptotic expansions with error control in Theorem 3. This allows us to 
prove that the integrability condition is not fulfilled. The proof of Theorem 4 therefore will 
be split into two parts: 

1. The first part consists in constructing a manifold M in the space of the parameters 
a, b, c, d such that if the eigenvalues for all Darboux points are real, then the param- 
eters belong to M. Then we produce a decomposition M = Mi U . . . Mk and study 
each manifold separately. For some of them, the corresponding potentials possess 
sufficiently many Darboux points to give a strong enough condition for integrability 
only using the Morales-Ramis criterion at order 1 (there could exist some resistant 
cases for which a higher variational equation is needed but without the phenomenon 
of arbitrary high eigenvalues like in 8 ). But for specific cases, this phenomenon occurs. 
It has already been noticed by Maciejweski in 9 who lets this specific case open. 

2. The second part will be devoted to these specific manifolds Mj where the Morales- 
Ramis criterion at order 1 is almost powerless. We use Theorems 2 and 3 to solve 
these hard cases. 

In 10 , the authors deal with a similar difficulty with the spring pendulum for which there 
is a discrete infinite set of parameters for which there are no obstructions to integrability 
at order 2. They also study third order variational equations, but then use analytic tools 
to study a sequence of monodromy elements, and finally prove that this sequence never 
vanishes. Thanks to our explicit expression via P-finite recurrences, such a problem can be 
analysed more systematically here. 

III. EIGENVALUE BOUNDING 
Definition 4. We will denote 



M 




r U (9) with U meromorphic and lit -periodic 
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Let V G M.. We denote by d(V) the set of Darboux points cofV with multiplicator — 1 and 
c 3 7^ 0. For c G d(V) we have Sp (V 2 V(c)) = {2, A} and we denote 



A(c) 



-oo otherwise 



Definition 5. We consider a subset E C M. and define 

A(E) = sup inf A(c). 

We sa?/ i/iai the problem of finding all meromorphically integrable potentials in E is a bounded 
eigenvalue problem if A(E) < oo. 

Remark 1. We have A(Ai) = oo because of the following family 

V(r, 8) = r- 1 ((1 + a) - 2ae w + ae 2id ) , aGl, 

for which only one Darboux point c = (1, 0) exists; the corresponding eigenvalue is A = 2a— 1. 
This proves that the family of potentials considered in Theorem 4 is an unbounded eigenvalue 
problem. 

Lemma 1. For a potential V G M. the Darboux points c such that C3 7^ can be written as 
c = (ci, C2) = (r cos(6 l o), r sin(6 l o)) with 8 being a critical point of U. The Darboux point c 
is non-degenerate if and only if U(6q) 7^ and in this case, the eigenvalues of the Hessian 
of V , evaluated at c, are 



if we choose the multiplicator of c to be —1. 

Proof. For V = r~ 1 U(8) the conditions (3) that c is a Darboux point are: 

r - 3 (- Cl U{6 ) - c 2 U'(6 Q )) = aci, 
r 3 (-c 2 U(6 ) + c 1 U'(6 ))=ac 2 . 

Assuming C3 7^ 0, it follows that U(8q) = -arg and U'(8o) = 0, which means that 80 is a 
critical point of U. Since c 3 7^ implies that r 7^ 0, we see that the case a = (degenerate 
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Darboux point) is equivalent to U(9q) = 0. Setting a = — 1 and U'(9q) = we get the 
Hessian matrix 

f(2c 2 1 -ci)u(e ) + c 2 2 u"(e ) Cl c 2 (3U(e )-u"(e )) \ 

„2 „2\TT(Q \ i „2tT"I 



V 2 V(c) 



'° \ Cl c 2 (3U(9 )-U"(9 )) (2cl-c 2 )U(9 )+clU"(9 )) 
whose eigenvalues are exactly those claimed above (using U{9q) = Tq). □ 

Recall that the potentials given by (6) are V(r,9) = r^ 1 U{9) with U(9) = a + be td + 
ce 2ld + de 3ld . We now assume that V possesses at least one non-degenerate Darboux point c 
with C3 7^ 0. After rotation, we can always assume that c = (1, 0) is a Darboux point. As 
shown in Lemma 1, it corresponds to a critical point for 9 = 0. Moreover, because this 
Darboux point is non-degenerate, we know that U(0) 7^ 0. Then by dilatation, we can also 
suppose that U(0) = 1 and get the following equations 

£/(0) = a + b + c + d=l, 
U'(0) = i(b + 2c + 3d) = 0. 

Solving these equations for c and d, yields the expression 

V a<b = r- 1 (a + be w + (3 - 3a - 2b)e 2l() + (2a + b - 2)e m 

for the potentials where a, b G C. 

Theorem 5. If V a ^ is meromorphically integrable, then it belongs to one of the following 
families 

E x = r~ l (-ho + 1 + be ie - be 2ie + \be m \ , 

E 2 = r- 1 (~k{k + l)e 3ie + h(k + l)e 2W - ^k 2 - ^k + 1 

E 3 = r- 1 (-h{k + l)e 2 ^ + \k{k + l)e w - \k 2 - \k + l) , 
_ 1 ( (s-6A 2 )A 2 m _ (3Ai + s - 3A 2 )A 2 2i9 



18(Ai + A 2 ) 6(Ax + A 2 

(6A 1 + s)A 2eie , -9AiA 2 - A 2 s + 18A X + 18A 2 - 3A 



2 s 



6(Ax + A 2 ) 18(Ai+A 2 ) 
where b G C and fcsN. T/ie quantities arising in E^ are 

s 2 = 6\j\ 2 + 6A1A 2 - 36AiA 2 , 

A*= i(A^-l)(Ar i + 2) + l (i = l,2), 

uitfi fci G N \ {0, 3} and A; 2 G N*. 
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Proof. For all non-degenerate Darboux points c = (7 008(6*0), 7 sin(6 l o)) the corresponding 
eigenvalue A satisfies 



(note that this condition is also satisfied if c is degenerate). We write U{6) = h a>b (exp(i9)), 
U'(6) = izh' ajb (exv(i9)) : and U"{9) = h a , b (exp(i6)) with 

h a:b (z) = a + bz + (3-3a- 2b)z 2 + (2a + b - 2)z 3 , 
~h ayb (z) = -bz - 4(3 - 3a - 2b)z 2 - 9(2a + b - 2)z 3 . 

So to find the eigenvalues of all Darboux points, one just needs to compute the following 
resultant which corresponds to the conditions (8): 

P a ,b(X) = res z (h a , b (z) - (A + l)h a ,b(z), h' ab (z)^ 




All the roots of -P a ,&(A) correspond to an eigenvalue of some Darboux point, except possi- 
bly in those cases (a, b) where P a ^ b vanishes as a polynomial in A or in the case where h' a b (z) 
has the root 0. 

Let us begin with the special cases. We compute the points (a, b) e C 2 for which P ajb = 
in C[A]. We find that it is the zero set of the ideal (2a + b — 2) n (a, 6). Moreover, the 
polynomial h' ab (z) has a zero root if and only if b = 0. So, all the specific cases belong to 
the zero set of (2a + b — 2) fl (b). First, for b = we find 



U"(6 ) - (A + l)U(9 ) = and U'%) = 



(8) 




(2a + b - 2)(6a + 2b - 6 + (A + 1))(-I8a6 2 - 66 3 + 186 2 + (A + 1) 
x (108a 3 + 108a 2 6 - 216a 2 + 36a& 2 + 108a - 108a& - 96 2 + 46 3 )) 



Qi = res z (h afl (z) - (A + l)h afi (z), h' afl (z)/z,z) 
= 216(a- l) 3 (6a-6 + (A + l)), 



and second, for b = 2 — 2a we get 



Q 2 = res z (K,2-2a(z) - (A + l)h a ^-2a(z), K ^ 2a (z), z) 
= -4(a-l) 2 (2a-2 + (A + l)). 



As we know that the eigenvalues should be of the form |(A; — 1)(A; + 2), k G N, we obtain 
the potentials E 2 and ^3 from these two cases. 
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Now for the generic case, we express a and b depending on the roots of P a ,&(A) and obtain 
the expression P 4 . Since it is not valid for k\ = k 2 = 0, we study this case separately and 
find the condition a = — 16 + 1, which gives E\. Note that fixing Ax = in E4 yields the 
potential E 2 , whereas Ai = 6 results in P3. The case k 2 = produces V = r _1 which already 
belongs to Pi. □ 

Corollary 1. With the same notation as in Theorem 5, we have A(Pi) = —1 and A(E 2 ) = 
A(P 3 ) = A(P 4 ) = 00. 

Remark 2. The types of E 2 , P3 and P4 differ fundamentally although they are all unbounded 
eigenvalue problems. This is because the dimension of P4 is 2 and the dimension of E 2 and 
P3 is only 1. Because of that, we could call P4 a doubly unbounded eigenvalue problem 
because it possesses two Darboux points whose eigenvalues can be independently arbitrarily 
high. Because of that, we will need to apply a third order integrability criterion simultaneously 
at the two Darboux points. The potential E\ has only one Darboux point with eigenvalue —1. 
This eigenvalue belongs to the Morales-Ramis table and so higher variational methods will 
be required, but only for this fixed eigenvalue (which is much easier). 

In the parameter space, we get 4 algebraic manifolds. For E 2 , P3, and P4, a tedious 
treatment with higher variational equations is required. For E\ we will be able to check 
integrability easily with Theorem 2. A similar procedure could be applied to any set of 
homogeneous potentials depending rationally on some parameters. Here computing power 
is the main limitation; in particular, because for typical problems, the number of parameters 
is much smaller than the number of roots which requires resultant computations and prime 
ideal decompositions. One should note that we have deliberately chosen a set of potentials (6) 
which is particularly difficult to treat. For most common problems (outside the general 
complete classification), these unbounded eigenvalue manifolds have small dimension (1 in 
the case found by 9 ) or even inexistent like in 8 or 11 . 

IV. HIGHER ORDER VARIATIONAL METHODS 

We will first recall some properties of the solutions of the first order variational equations. 
After diagonalisation and in the integrable case, the equation is the following (after fixing 
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the energy E = 1) 

2t 2 (l + t)y"(t) - ty'(t) - \{n - l)(n + 2)y(t) =0 (n G N). (9) 

After the change of variables t — > {t 2 — this equation becomes 

(i 2 - l)y"{t) +4ty'(t) - (n- l)(n + 2)j/(t) = (n e N). (10) 

A basis of solutions is given by (P n , Q n ) where P n are polynomials in t (for n > 1) and the 
functions Q n are 



Q„(f) = Pn(t) J 



(t 2 -l) 2 P n (t) 2 ^ 

The functions Q n are multivalued except for n = which will be a special case. Indeed, the 
Galois group of (9) in this case is Id instead of C. 

The polynomials P n can be computed using the Rodrigues type formula 

1 an-l 

Pn(t) = W^ldt^^ - ^ (^!) ( n ) 

which also gives a normalisation for their leading coefficient. The functions Q n can be 
written as 

1\ , W n (t) 

where W n are polynomials given by 

1 +2 



Q n (t) = e n P n {t) arctanh ( 7 ) + -^f (n > 1) (12) 



^) ^- i ^- i M-^^: M ^- t - + 



r(* + ±)' 



W 2 fc+l(t) 



2fct(2fc + 1) arctanh(t) aFWl - k,k+ §, §,t 2 
(fc!) 2 

(_l)fc( t 2 _ !) + i)(2A; + l) 2 Fi(§ -k,k + 2, It 2 ) 



2 4fc + 2 V r(fc + |) 2 

3 1 
2' 2' 



2arctanh(t) 2 Fi(-k, k + |, §, t 2 



(fc!) 2 ; 

and e„ is a real sequence given by 

n(n + 1) 
6n = 4™(n!) 2 ' 

Conventionally, we will take for n = 0: 

m = ^y, Qo(f) = ^ =rr 
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Lemma 2. The functions P n (t) and ^Q n (t) satisfy the differential equation (10) and the 
three-term recurrence 

An{n + l)(n + 2)y n (t) - 2t(n + 2)(2n + 3)y n+ i(t) + (n + 3)y n+2 (t) = 0. 

Proof. Given the explicit expressions (11) and (12) we can use holonomic closure properties 
to derive the differential equation resp. recurrence they satisfy. We first express (11) as 

, * (n-l)\ r(u 2 -l) n , 



by Cauchy's differentiation formula. By the method of creative telescoping we obtain the 
differential equation and the recurrence (this calculation was carried out by the software 
package HolonomicFunctions 12,13 ). Similarly we can apply holonomic closure properties to 
the closed form expression (12). □ 

Lemma 3. (proved in 14 ) Let F G C(zi) [z 2 ] and f(t) = F(t, arctanhf|^ . We consider the 
field extension 

K = c(t, arctanh Q j , J f dtj 
and the monodromy group G = a(K,C(t)). If G is abelian, then 

— Res t=00 F (t, arctanh (-^\ + a^j =0. 

Proof. We will consider two paths, the "eight" path o~i around the singularities —1 and 1, 
and the path cr 2 around infinity. At infinity, the function F(t, arctanh (j) + a J will have a 
series expansion of the kind 

J F (t, arctanh (^j + cxj dt — ^ a n (a)t n + r(a) hit 

because the function arctanhQ) has a regular point at infinity. Let us now consider the 
monodromy commutator 

-i -JL. JL. 

We have that crp r (/) = F(t, arctanh (j) + and cr 2 (lnt) = lnt + 2in. We deduce that 

ff(/) = / + r(/9)-r(0). 

This r(/3) corresponds to the residue of F(t, arctanhQ^ + f3^j at infinity. If the monodromy 
is commutative, then the commutator a should act trivially on /. This is the case only if 
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r(/3) — r(0) = for all G 2?7rZ. The function r is a polynomial in (5, so r((3) — r(0) = for 
all (3 e C. From this the claim follows. □ 



In the following, we will also need to use the next lemma which is a kind of reciprocal 
version of Lemma 3. 



Lemma 4. (proved in r ) We consider 



3 /i 
F(t) = Hi(t) arctanh 

i=0 



t, arctanh ( | 



w'i/i i^o? • • • , i?3 G C[t] . // i/ie conditions of Lemma 3 are satisfied, then 

• If Res t=oc F(t) = 0, then f F dt e C 

• // Res t=oc ^ 0, then JFdteCk, arctanh (|) , In {f - 1) 

Theorem 6. Let V be a homogeneous potential of degree —1 in the plane. We suppose 
that c = (1, 0) is a Darboux point of V with multiplicator —1. If the variational equation is 
integrable at order 2 then 

Sp(v 2 \/( C )) = {2,i(p-l)(p + 2)}, peN, 

and for odd p we have = 0. 

This theorem is in fact a particular case of Theorem 2 in 14 for which the three indices 
i,j,k are equal. 

Remark 3. Because the constraint appears only for odd p, the variational equations of 
order 2 give no constraint for even p. Hence this is not sufficient for proving non-integrability 
for an unbounded manifold. 

14 



Integrability conditions for homogeneous potentials 



V. PROOF OF THEOREM 2 

Proof. The variational equation at order 3 is given by 

X - 2 X + 1 a Y Ab 

A a b c 

X2 = J^ 2 + J^ 1 + J^ 1 + ^ Z 

= 2Y 1>2 
Yl ' 2 = ^ Yl ' 1 + ^ Z3 + Yl ' 3 

Y 2y \ = *2,2 + *2,3 



3 



Y 2j2 = 


03^ 


46 , 

,Z 3 + Y 2 ± 

30 5 ' 


^2,3 — 


+ Y 2A 


^2,4 — 




-^Z + ±Y 22 
30 5 3 ' 


z = 







where A = \{n — l){n + 2). The coefficients a, b, c correspond to the following derivatives 

" = <W< C >' b -2M v{c) ' C =SW? [C) < 

and the others are given using the Euler relation for homogeneous functions. A complete 
procedure to build these equations is given by 15 . The functions Y^x and Y 2 ^ are solutions of 
a system of linear differential equations with an inhomogeneous term, and the homogeneous 
part is in fact a symmetric product of the first order variational equation. Here, we already 
put to zero terms that we think in advance they will not produce integrability constraints. 
As before, we use the change of variables <fi(t) — > (t 2 — l) -1 . 

We choose Z(t) = Q n and compute the solution for X 2 of the above system. We first 
remark that X 2 is in the Picard-Vessiot field, so it is also the case for its derivative. We now 
perform integration by parts and see that one term is already in the Picard-Vessiot field, 
and the other is 

J 2(t 2 - l) 2 (a 2 tP n Q n h + Ab 2 P 2 Q n I 2 + c(t 2 - l)Q A n ) dt (13) 
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where 



' / ^ t 2 (t 2 - l) 2 + (t 2 -l) 2 P n 2 y | 

/ ((t 2 - l) 2 g3 + _ S _2 J3 _ I S (f _ l)4p n Q 2 (P n g n - g n p n )dt) dt 
h ~ 1 " (t 2 - l) 2 P n 2 dt 



4 = / t(t 2 - i) 4 g 2 (p n g n - Q n P n ) dt 



Let us now study this expression term by term. We begin with the third summand of (13) 
which is 

2c / (t 2 -l) 3 Q*dt. 



It has already the form of Lemma 3. So as in the proof of Lemma 3, the monodromy 
commutator will be computed using 

Res t=00 (t 2 -l) 3 (g n + e n «P n ) 4 . 

Now look at the term in b 2 . It is not as complicated as we could think because of the 
following relation 

PnQn - PnQn = " 1)~ 2 Vn G N 

which is linked to the Wronskian of Equation (10). Thanks to that, the term in b 2 can be 
written as 

and then using integration by parts, this gives 



IQb 2 J Qlit 2 - I) 2 J P n Q 2 n (t 2 - ifdtdt - 8b 2 J P n Q 2 n (t 2 - l) 2 dt J Qlit 2 - ifdt 
Now by Lemma 4 we have for all even integers n > 1: 



J P n Q 2 n (t 2 - Ifdt, J Qlit 2 - Ifdt G C(t) 



arctanh 



So we are integrating a polynomial in arctanh with rational coefficients, and this corresponds 
to the hypotheses of Lemma 3. The second term does not provide any monodromy, so we 
only have to study the first term and thus the sequence 

Res t=00 {Q n + e n aP n f(t 2 - \f J P n (Q n + e n aP n f{t 2 - if. 
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Now we look at the term in a 2 . It can be simplified to 



r {t 2 -i) 2 Qlt , J(t 2 -i) 2 Qltdt j + r f(t 2 -i) 2 Q 2 n tdt , 

We now use again integrations by parts (recall that P 2 = At): 

8a 2 J(t 2 - l) 2 Q 2 n Q 2 J(t 2 - ifQltdt - 8a 2 J(t 2 - l) 2 Q 2 n t J(t 2 - l) 2 Q 2 n Q 2 dt. 



To conclude we can again use Lemmas 3 and 4. We first prove that 



\/n £ 1 J P 2 Q 2 n (t 2 - l) 2 dt, J Q 2 n Q 2 (t 2 - l) 2 dt G C(t) 



arctanh 



The case n — 1 corresponds to A = 0, for which we have always the coefficient a = 0. Now 
we make a final integration by parts which gives 

16a 2 J{t 2 - l) 2 Q 2 n Q 2 J{t 2 - IfQltdt dt - Sa 2 J{t 2 - l) 2 Q 2 n tdt J[t 2 - l) 2 Q 2 n Q 2 dt. 

Thanks to that, we get a constraint of the form given by Theorem 2 and the coefficients are 
given by (multiplying them by e~ 2 for further simplifications) 



h(n) = (a 3 ) 2e~ 2 Res t=oc \^t 2 - l) 2 (Q n + e n aP n ) 2 {Q 2 + e 2 aP 2 

x J(t 2 -l) 2 {Q n + e n aP n ) 2 P 2 dt^, 
f 2 (n) = (a 3 ) 2e~ 2 Res t=OQ ((t 2 - l) 2 (Q n + e n aP n ) 3 

x J(t 2 -l) 2 (Q n + e n aP n ) 2 P n dty 
h{n) = (a 3 ) U- 2 Res t=00 ((t 2 - l) 3 (Q n + t n aP n ) A ) , 



(14) 



(15) 



(16) 



where (•) denotes coefficient extraction. In fact, only the coefficient of a 3 appears in these 
residues. We need not to prove this fact, because we simply select the coefficient of a 3 , 
ignoring the question whether the other coefficients are zero or not. 

We now look at the case n — 0. All our previous calculations are also valid in this case 
except those involving Lemma 4 because we only have 

r 



J P Q Q 2 (t 2 - l) 2 dt, J Q 3 (t 2 - l) 2 dt g C(t) 

J P 2 Q 2 (t 2 - l) 2 dt, J Q 2 Q 2 (t 2 - l) 2 dt G C(t) 



arctanh [ ™ ) , ln(t 2 — 1] 



arctanh 
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So, the coefficients in a 2 , c are also 

2Res t=00 (V - l) 2 Q 2 Q 2 J (t 2 - l) 2 P 2 Q 2 dtj , 

gRes 4=00 ((i 2 -l) 3 Q*). 

We find that these residues are both 0, and so the corresponding integral does not provide 
any additional monodromy. The case of the coefficient in b 2 is a little more difficult because 
the integral does not satisfy the conditions of Lemma 3. After an explicit computation, we 
arrive at the following integral 

J -^—^ (-t arctanh (^j - ^ In (t 2 - i^J dt = 
^ln(2)ln(t-l) + ^dilog(t + l) + 
gln(t + l) 2 + ^ln(t + l)ln(t-l)-gln(t-l) 2 . 

All the terms are in C[t, arctanh yjj ,ln (t 2 — 1)] except one, namely the dilogarithmic term 

d 1 log(t + l) = / ln(t + 1) rft. 

With the same idea as in Lemma 3, we see that this term has a noncommutative monodromy 
because of the following residue in 

ln(t + l) + a 
Res t=0 = a 

which depends explicitly on a. So, for n — 0, the integrability condition at order 3 is in fact 
just b 2 = 0. 

□ 



VI. HOLONOMICITY AND ASYMPTOTICS 

In this section we are going to derive P-finite recurrences (i.e., linear recurrences with 
polynomial coefficients) for the sequences fi(n), f2(n), and fzin) that appeared in section 
V. The methods that we employ are based on Zeilberger's holonomic systems approach 16 . 
The recurrences presented below were computed with the method of creative telescoping, to 
which a brief introduction is given below (see 12 for more details). 
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Let S n denote the forward shift operator in n, i.e., S n f(n) = f(n + 1), and D x the deriva- 
tive w.r.t. x, i.e., D x f(x) = f'{x). The method works for the class of holonomic functions, 
which in short are (multivariate) functions that are solutions of maximally overdetermined 
systems of linear difference and differential equations with polynomial coefficients. The set 
of all equations which a given holonomic function satisfies forms a left ideal (we call it 
annihilating ideal) in some Ore algebra of the form 

C(m,n, ...,x,y,.. 

The nice fact about holonomic functions is that this class is closed under certain operations 
(addition, multiplication, certain substitutions, definite summation and integration) which 
can be executed algorithmically: given the defining systems of equations for two holonomic 
functions / and g, there are algorithms to compute a holonomic system for / + g, f ■ g, etc. 

For computing integrals (or residues), the method of creative telescoping makes use of the 
fundamental theorem of calculus. Consider a definite integral of the form J% f dx where the 
integrand / depends also on some other (discrete and/or continuous) parameters. We need / 
to be holonomic, i.e., there is some left ideal I of annihilating operators in the corresponding 
Ore algebra O. The idea is now to come up with an operator A + D X B e / where A, B G O 
and A does not depend on x and D x (the concept of Grobner bases 17 plays a crucial role in 
this step). Then after integration we get, 



P 



fdx + 

J a 



i 6 



Qf 



o. 



in other words, we found a (possibly inhomogeneous) equation for the integral in question. 
The examples below will demonstrate this methodology clearly; we start with the simplest 
one, the sequence fziji). 

Lemma 5. The sequence fz(n) given in (16), satisfies the P-finite recurrence 
(An + 11) (4n + 9)(n + l) 3 (n + 3) 2 f 3 (n + 2) - 

(2n + 3)(16n 6 + 144n 5 + 515n 4 + 930n 3 + 888n 2 + 423n + 81)/ 3 (n + 1) + 
(An + 3)(4n + l)(n + 2fn 2 f 3 (n) = 0. 



subject to the initial conditions 



M 1 )- ~105' /3(2) -"385' 
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Proof. It is an easy exercise to compute the first values of fsfo) explicitly with a computer 
algebra system. Thus we basically have to derive the recurrence. For this purpose, we 
compute an annihilating ideal / for (t 2 — l) 3 (Q n + e n aP n ) A which is the expression in the 
residue (16). For this purpose we apply holonomic closure properties (note that Q n + e n aP n 
satisfies the same equations as Q n itself). The resulting Grobner basis is too large to be 
printed here, namely a full page of equations approximately. It is represented in the Ore 
algebra C(n, t)(S n , D t ). In the next step we make use of a special algorithm 18 for computing 
a creative telescoping operator 

A{n,S n ) + D t B{n,t,S n ,D t ) el 

(its existence is guaranteed by the theory of holonomy). Because we are dealing with a 
residue we can forget about the part B and find that A annihilates the residue. In order 
to obtain fz(n) we need to multiply the residue with 2e~ 2 , which can be done again by 
closure properties. The resulting operator represents exactly the above recurrence. All these 
computations were done with the above mentioned package HolonomicFunctions 12,13 . □ 

Lemma 6. The sequence f\(n) given in (14) satisfies the P-finite recurrence 

(An + ll)(4n + 9)(n + 4) 2 (n + l) 3 (4n 2 + 8n - 9)/i(n + 2) - 
(2n + 3)(64n 8 + 768n 7 + 3580n 6 + 8028n 5 + 8113n 4 + 

834n 3 - 4863n 2 - 3276n - 648)/i(n + 1) + 
(An + 3)(4n + l)(n + 2f(n - l) 2 (An 2 + 16n + 3)/ x (n) = 

subject to the initial conditions 

/l(2)= 1155' /l(3) = 2145- 

Proof. The proof is based on the same ideas as in Lemma 5, except that the expression of 
which we have to take the residue is more complicated. In particular, an indefinite integral 
occurs (recall that indefinite integration is not among the holonomic closure properties) and 
it is not clear a priori how to choose the integration constant such that the result is again 
holonomic. We start by computing an annihilating ideal I for 

F(n,t) = (t 2 -l) 2 (Q n + e n aP n ) 2 . 
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Thus for all A G / the operator AD t annihilates the indefinite integral J F(n, t) dt. Addition- 
ally, from a creative telescoping operator A + D t B G / we can derive more such annihilating 
operators. Let J denote the annihilating ideal for B(F) which can be obtained by holo- 
nomic closure properties. Then for every C G J, the operator CA annihilates the indefinite 
integral as well. Altogether we obtain a zero-dimensional annihilating ideal for / F(n, t) dt, 
and continue as in Lemma 5. □ 

These recurrences in Lemmas 5 and 6 are irreducible (in the sense that the corresponding 

operator cannot be factorized), and so we are not able to find closed forms for f\ and f 3 . 

The recurrence for /2(2n) is given by a third-order recurrence with polynomial coefficients 

of degree larger than 50, which we do not state here explicitly. The initial conditions are 

, / x 16 , / iX 184 , . . 38308 
/a(2) = 777^ /2(4) = 7777777"; / 2 (6) 



1155 J y ' 183141' J y ' 181081875 

This recurrence is reducible and possesses a hypergeometric solution 

8vr 2 r(n + l)r(5/6 + n) 2 T(l/6 + n) 2 Y{nf 
^ >25T{n + 2/3) 2 r(3/2 + nfT(l/2 + n)r(4/3 + n) 2 

but because ^2(2) 7^ 0, the recurrence for f 2(^11) cannot be reduced. 

We are interested in a practical way to apply the third-order variational equation. To 

do this, these recurrences are not enough, since we need closed forms. As these closed forms 

do not exist, we will instead produce closed form expressions which approach f ly f 2 , and f 3 

with a controlled relative error. In the following, we will denote the harmonic numbers by 

n— 1 1 

H{n) = £ - 

i=l 

Definition 6. Let us consider an operator L G C(n,S n ), in other words L represents a 
linear recurrence with polynomial coefficients. We will say that L is regular at infinity if for 
all solutions u (i.e., Lu = 0) there exist a G Z, /3 G N, and 7 G C such that 

u(n) ~ ^n a H(nf for n — > 00. 

Theorem 7. Consider L G C(n,S n ) of order k and assume that it is regular at infinity. 
Then for all p G N and for all u solution of Lu = 0, there exists a function F G C(n)[H(n)] 
with degree in H{n) less than k — 1 such that 

/ H(n) k - 1 \ 



u(n) = Fin) + O 



n p 
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This theorem is directly implied by the theorem of Birkoff given in 19 , which gives a 
general form of an asymptotic expansion which is always possible. In our case, we will only 
use what we call the regular case, which in a Birkoff expansion corresponds to not having 
an exponential part. 

Definition 7. Consider a function f : N — > M. and a function F 6 M.(n)[H(n)]. We say 
that F is an approximation of f with relative error e at rank Uq if 

f(n) 



F(n) 



- 1 



< e Wn > uq. 



We consider p functions fi, ■ ■ ■ , f p '■ N — > K and approximations Fi, 
with relative error e at rank n . We define the error amplification factor A by 

v 

A = min Mell such that 



F p G R(n)[H(n)] 



E fi(n) 
i=i 



£ Fi{n) 

i=l 



- 1 



< Ae Vn > no 



> . 



Lemma 7. We consider p functions fi, ■ ■ ■ , f p '■ N — > R and approximations . . . , F p 6 
M.(n)[H(n)] with relative error e < 1 at rank n and A their amplification factor. Then 



E \Fi(n 
A < max % ~~ 

n>no 



E Fi(n) 
i=i 



Proof. The lemma is equivalent to prove that 



E fi(n) 
i=i 



E Fin) 



i=i 



< emax 

n>no 



E \F(n) 

i=l 



E Fin) 



i=i 



So one just needs to maximize the left hand side. We already know that \ f%{n) / F^n) — 1 1 < e. 
So depending on the sign of fi(n) we replace fi(n) by (1 — e)Fj(n) or (1 + e)i^(n). We then 
expand 



E /i(n) 

i=l 

E Fin) 
i=i 



- 1 



< 



e £ sign(/i(n))Fi(n) 

i=l 

E FiCn) 
i=i 



< 



e E 

i=l 

E F(n) 
i=i 



< e max 

n>no 



E \F(n) 

i=l 



E Fi(n) 



i=i 



using the fact that fi(n) and Fi(n) have always the same sign for n> n (because e < 1). □ 
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In practice, we first check that the sign of the functions Fi(ri) and their sum does not 
change for n > n and then we prove a majoration of the resulting expression in M.(n, H{n)). 
So all comes down to prove that some polynomial in IR[n, H(n)} does not vanish for n > n . 
This can be done by first making an encadrement of the function H(n) and then prove that 
the corresponding bivariate polynomial does not vanish on a particular algebraic subset. 
Such a problem can be algorithmically decided. 

Theorem 8. Consider the recurrence equation 

u{n + 1) = A(n)u(n) Wn G N, A(n) G M P (C) (17) 

Consider \\-\\ a matricial norm and R(n) the resolvant matrix of equation (17). Assume that 

oo 

m(oo) = J2\\Mj)-i p \\ <i- 

j=0 



Then 



Proof. We write 



Let us pose 



n—1 n— 1 

R(n)=l[A(i) = l[((A(t)-I p ) + I p ). 

i=0 i=0 

M(n) = J2\\A(j)-I p \\. 

3=0 

We want to prove a majoration of the type 

\\R(n) - I p \\ <CM(n) (18) 

with a suitable constant C > 0. For n — 1, this is true with (7 = 1. Let us prove 
equation (18) by recurrence: 

m = nV« - j P ) + j P ) = w - 1) - j P ) ff^w + ff^)> 

i=Q i=0 i=0 

i?(j) - R(j - 1) = - 1) - I p ) = - 1) - - 1). 
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Then we sum these equations for 1 < j < n which produces 



\\R(n)-I p 



n-1 

3=0 

n-1 

< J2\\ A U) - i P \\ \\RU)-i P \\ + \\m-ip 

j=0 
n-1 

< M(n) + J2\\ A (j) - I P \\ CM{j) = M{n) + CM(n) 



j=0 



< (1 + CM{oo))M{n) 



using the fact that M(n) is a growing sequence. So the recurrence property is proved if 
C < 1 + CM(oo) which is equivalent to C > (1 — M(oo)) -1 > 1. So this proves that 

M(n) M(oo) 



||i2(n) 



',11 < r 



< 



M(oo) ~ 1 - M(oo) 



which proves the theorem. 



□ 



The main application of this theorem is to compute a sequence with controlled error. Let 
us take an operator L G lR(n, S n ) regular at infinity. We can then compute an asymptotic 
expansion of the resolvant matrix of L, and an error matrix which will satisfy an equation 
like (17). Then for an n e N, we can apply Theorem 8 for the shifted sequence u(n + no), 
and the majoration M(oo) will become very small for no big enough, giving us that the error 
is always lower than some explicit bound. This has very important consequences for the 
application of the higher variational method. In particular, it becomes possible to rigorously 
prove that a sequence of potentials with the unbounded eigenvalue property does not satisfy 
integrability criteria for A large enough, and thus coming back to a bounded eigenvalue 
problem. 



VII. APPLICATION AT ORDER 2 

We now apply the second order criterion to our example. We begin with the case E4. 
Before we state the corresponding theorem, we need a preparatory lemma concerning the 
solutions of a certain Diophantine equation. 

Lemma 8. The set of solutions (/ci,fc 2 ) G N 2 of the Diophantine equation 

R(k u k 2 ) = k\k\ + k 2 k\ - 75A; 2 - 75^ + k 2 h - 27k 2 + k\k x - 27kj = 
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is given by {(0,0), (6,14)}. 

Proof. We begin by proving that for k 2 > 50, the condition R = implies 4 < k\ < 5, 
and similarly, for k\ > 50, we have 8 < k 2 < 9. These statements can be written as logical 
expressions involving polynomial inequalities 

\fh\fk 2 : (k x > A k 2 > 50 A R[k u k 2 ) = 0) =^ 4 < fc x < 5, (19) 
\fh\fk 2 : (h > 50 A k 2 > A i2(ifei, fc 2 ) = 0) 8 < k 2 < 9. (20) 

Such formulas can be proven routinely with quantifier elimination techniques like cylindrical 
algebraic decomposition 20 . Indeed, applying the Mathematica command CylindricalDe- 
composition to the above formulae reveals that they are true. Therefore, there are no 
integer solutions for k% > 50 or k 2 > 50 and an exhaustive search delivers exactly the 
solutions claimed above (Figure VII. 

However, if we want to prove (19) and (20) "by hand" (let's consider the first one for the 
moment), we have to look at the largest real root of the polynomial 

res fcl U(k t , k 2 ), 9 ^^- ) £(4, h) R(5, k 2 ). 
We find that this root is smaller than 50 (using real root isolation) and that the limit 

lim n(k 2 ) = -\ + 

k,2— »oo 2 2 

is between 4 and 5, where n(k 2 ) denotes the positive solution of R(ki, k 2 ) = regarded as 
an equation in k\. The implication (19) follows, and (20) can be proven analogously. □ 

Theorem 9. We consider the potential E± given in Theorem 5. If the variational equation 
near all Darboux points is integrable at order 2, then the corresponding eigenvalues are 
integers of the form X = (21 - 1)(Z + 1), I G N. 

Proof. We use the notation U = rE^ from Theorem 5. The condition U'(6) = yields the 
two Darboux points 

d : e" = 1, 

* ■ - ^ (21) 
s - 6A 2 

There are singular cases of the second equation, namely for s + 6Ai = or s — 6A 2 = 0. After 
solving and replacing, we find that these cases correspond exactly to ki = and k± = 3, 
which were excluded from £4. 
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We now compute the third derivative of V, evaluated at the two Darboux points c\ and c 2 
given by expression (21): 

d 3 V _ iAi(g + 15A! + 9A 2 ) 
dql (ClJ ~ Ai + A 2 

d 3 V _iA 2 (s-15A 2 -9Ai) 
M [C2)= 3(A 1 + A 2 ) ~- 

In the case (k\,k 2 ) both odd, both derivatives should vanish. We solve the system and we 

find Ai(k 2 + l)k 2 = 0. This is impossible for odd values. In the case k\ odd k 2 even, the first 

one should vanish, and in the case k\ even k 2 odd the second one should vanish. We get the 

equations 

klikx + \f{k\k\ + k\k x - TJk\ - 27k 2 - 75h + k 2 h - 75kf + k 2 kf) 

12(k 2 2 + k 2 + h + kl) 
k\{k 2 + \y{k\k\ + k\k 2 - 27fc? - 27fci - 75fc 2 + fcifc 2 - 75jfc| + fcifcl) 1 ' 

12(kf + h + k 2 + kl) 

These two conditions are symmetric. The first terms can never vanish because we have 
k\ odd for the first one and k 2 odd for the second one. To conclude, we need to look at the 
last term, which corresponds to a Diophantine equation, and to prove that this equation 
does not have a solution with k\ odd and k 2 even. 

With Lemma 8, we have no solutions from the second term where k\ and k 2 have different 
parity. We conclude that all the possibilities left are for k\,k 2 even. □ 

It is well known that Diophantine equations in general cannot be solved (Matiyasevich's 
theorem). This means that Lemma 8 is a lucky case, although not trivial to prove. We 
therefore should remark that the study of this equation is not absolutely mandatory. We 
could simply skip it, assume that it is satisfied and continue further to the third-order 
condition. This condition would add two additional equations in k\ and k 2 and thus would 
allow to solve the problem in all generality. 

Here we are in a special case. A Diophantine equation R(k\, k 2 ) = can be solved only 
using real algebraic geometry in one of the following cases: 

1. The set 7? _1 (0) D M +2 is compact. In this case we only have a finite number of points 
to test. 

2. The set -R _1 (0) HlR +2 is not compact but all infinite branches are asymptotes and the 
corresponding asymptotic straight lines have a rational slope. In this case, either R 
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FIG. 1. Graph of R 1 (0). The graph R 1 (0) n M + is not compact but the infinite branches are 



k2 + I - 3V301 = 0. 

is homogeneous and has an infinite number of solutions, or the integer solutions can 
be bounded: when approaching infinity, the infinite branch of _R _1 (0) comes closer to 
the asymptotic line without touching it; for rational slope, there is then a nonzero 
infimum for the distance between the asymptotic straight line and integer points). 

The first case can be considered to be part of the second one with no asymptotes at all. In 
Lemma 8, we encounter the second case. 

Remark 4. The potential corresponding to ki,k 2 = (6,14) is the following (with the good 
choice of valuation for the square root): 



This potential has two Darboux points, it is integrable at order 2 near these two Darboux 
points and we have also that the third derivative near one of the Darboux points is zero 
(which is not needed for integrability at order 2 but gives interesting properties in practice 
at order 3 ). 

Theorem 10. Among the potentials in the families Ei, E 2 , E 3 , if a potential V is meromor- 



asymptotes with rational or vertical slopes; here the asymptotes 
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phically integrable, then it is of the form (after multiplying by some constant factor): 

V= X - (-h(2k + l)e m + ^k(2k + l)e 2ie - ^(2k 2 + k - 6) 

V = 1 (- 1 k(2k + l)e 2W + k(2k + l)e ie - l (2k 2 + k - 2) ), 
r V 2 2 / 

for fceN. 

Proof. The potentials E 2 and E 3 possess only one Darboux point. The corresponding po- 
tentials are 

E 2 : V = r- 1 (~k(k + l)e™ + ^k(k + l)e 2 ' ld - ^k 2 - ^k + 1 
E 3 : V = r' 1 (-\k(k + l)e 2ld + h(k + l)e %e - h 2 - + l) . 

We know that if k is odd, we have an additional integrability condition at order 2. We find 
that 

d 3 V 5 d 3 V 3 

— (c) = -ik(k + 1) for £ 2 , —(c) = -zA;(A; + 1) for E 3 . 

oq 2 I Oq 2 ^ 

These terms should vanish. This is never fulfilled for odd k. The sequence of potentials given 
by Theorem 10 corresponds exactly to the cases of even k (for which there is no condition for 
integrability at order 2). At last, we have the potential E%. The corresponding eigenvalue is 
always —1, so it is always integrable at order 2. At order 3, we know that the integrability 
condition is U^ 3 \0) = 0. We get 

U {3) (0) = -2ib 

So the only possibility is b = and this corresponds to the potential V = r _1 . This potential 
is integrable and already belongs to the family described by Theorem 10. 

□ 

VIII. APPLICATION AT ORDER 3 

We will now prove Theorem 3, building an algorithm to prove it. 
Proof. The scheme of the proof is the following 

• First we prove that the recurrences for fi, f 2 , f 3 are regular at infinity. 

• We then produce a series expansion Ri(n) at infinity at an order high enough of the 
resolvant matrix Ri(n) associated to these recurrences. 
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We then write Ri(n) = Ri(n)Ri(no)~ 1 Ri(no)Ei(n) for a large enough fi GN and build 
a recurrence of the form (17) whose resolvant matrix is Ei(n) (after change of basis), 
which will be denoted by Ei(n + 1) = Ai(n)Ei(n). We have moreover that Ei(no) is 
the identity matrix. 

As Ri(n) is a good approximation of Ri(n) when n — > oo, the matrix A^n) will tend 
to the identity matrix when n — > oo. Using Theorem 8 with a shift in the indices, 
we will have that 



E 

/="0 

oo 

i - E \\Mj) - i\ 



\Ei(n) - I\\ < Vn>n 



j=n 

• If we have chosen an expansion order and uq large enough, this sum will be finite and 
small, and thus will give us an approximation of Ri(n) by Ri(n) with relative error 
control. The expressions in Theorem 3 follow. 

For fs(2n), we find the following asymptotic expansion (a high order makes up the com- 
putation easier for error control) 

(\ 1 25 35 183 \ 
Cl Vn* ~ n? + 32r^ ~ 64^ + 512n 8 / + 

° 2 VV16n 4 ~ 16n 5 + 512n 6 _ 1024n 7 + 8192n 8 / ^ + 
1 1 19951 7507 96541 58151 

~ 2r^ + 46848n 4 ~ 46848n 5 + 1499136n 6 ~ 2998272n 7 

This proves by the way that the recurrence for /3(2n) is regular. We do the same for /i(2n) 
and /2(2n) and we find that they are regular too. We then find a majoration of the norm of 
the error matrix A 3 (n) 

< 9975 29925 H (n) 9975 29925 H(n) 
" 3 ^"°° ~ 256n 6 + 4096 rfi~ + 256n 8 + 4096 n s 

We choose now no = 100. We majorate the sum of this majoration beginning at n — 100. 
We find a majoration of this sum by 



E WMn) Hoc < 4.84522 x 10" 

n=100 



9 



4.84522 x IP' 9 
1 v ; " ~ 1 -4.84522 x 10~ 9 
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(an explicit rational number). We then compute the recurrence up to n = 100, and then 

produce an encadrement (with error less than 10~ 10 ) of the result with rational numbers. 

Although it is not mandatory in theory, in practice recurrences tend to produce very large 

rational numbers, whose size grows linearly with n, and thus are impractical to manipulate. 

This gives us the coefficients ci,c 2 with a good error control: 

883919839 1740684681 



c l = ~ 77777777777777777777 > c 2 



274877906944' 8589934592 
We then compute the error amplification of the sum, and find that it is less than 33/32. As 
the resulting expression is too complicated to manipulate for applications, we only keep the 
terms up to order 3 and prove that this new approximation has a relative error less than 
10~ 5 . The expressions for f\ and f 2 are found with a similar way, with the exception that 
at the end, to produce a sufficiently simple and accurate formula, it is not sufficient to keep 
the terms up to order 2 (after there is a H(n) that we want to avoid), so we need to add 
a term of order 3 (without H(n)) with a well chosen coefficient such that the error stays 
below 10~ 5 (else the result is only accurate to 10~ 3 ). □ 

Theorem 11. The third order integrability conditions for the families 

V = - (-\k(2k + l)e 3W + h(2k + l)e 2W - \{2k 2 + k - 6)) 

V=- (~\k(2k + l)e 2W + k(2k + l)e w - ^(2k 2 + k - 2)) 
r \ z a j 

where k G N* , are 

9(k + l) 2 {2k - l) 2 f\(2k) = 25k 2 {2k + l) 2 / 2 (2fc) + (66A; 2 + 33A; - 9)/ 3 (2ife), 
9(k + l) 2 (2k - l) 2 fi(2k) = 9k 2 (2k + l) 2 f 2 (2k) + (42A; 2 + 21k - 9)f 3 (2k), 

respectively. They are never satisfied. 

Proof. We replace fi(2k),f 2 (2k),f 3 (2k) by their approximations, and then compute the 
error amplification. It is less than 33/32, and the resulting expression does not vanish for 
k > 100. For k < 100, we make exhaustive testing and we do not find any solutions. For 
the second equation, we do not find any solution either. □ 



Theorem 12. We consider the family of potentials E4 

(s-6A 2 )A 2 3 -g _ (3_A 
18(A! + A 2 ) 6(Ax + A 2 ) 



£ 4 . v = r -i( (s-6A 2 )A 2 c3 .g _ (3A! + s-3A 2 )A 2 c2 . e | 



(6A1 + s)A 2 ie -9AiA 2 - A 2 s + I8A1 + 18A 2 - 3A 



2 s 



6(A! + A 2 ) 18(Ai + A 2 ) 
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with 




The third order integrability condition for E4 is of the form 



Qk lM (fi(h),f2(k 1 ),f 3 (k 1 ))=0 



Qk 2 , kl (fM,f 2 (k 2 )J 3 {k 2 ))^0 



where Q is a quadratic form depending polynomially on k\ and k 2 . 

Proof. We use Theorem 2 and compute the derivatives of the potentials in the family 
These derivatives depend rationally on k\, k 2 , and s. As there are two Darboux points, we get 
two conditions (d), (C 2 ) linearly dependent on fi(h), f 2 (h), f 3 (h) or fi(k 2 ), f 2 (k 2 ), f 3 (k 2 ) 
respectively for each Darboux point. To remove the quadratic extension s, we make the 
product (Ci) x subs(s = — s, (Ci)) and (C 2 ) x subs(s = — s, {C 2 )). The fact that in the 
potentials of £4, the two parameters Ai and A2 play a symmetric role produces the two 
conditions Qk lt k 2 — an d Qk 2 M = 0- ^ 

Remark 5. The conditions Qk 1 ,k 2 ,Qk 2 ,k 1 ar & n °t equivalent to the conditions (Ci),(C 2 ). 
We can solve (C\) in the quadratic extension and get for example that s should be rational 
because fi, f 2 , fs are always rational (this can be proven even without the P-finite recurrences 
since they correspond to a particular term in the series expansion of rational expressions in 
t,P n (t),Q n (t)). We get that 



if some generic condition depending on the fi(k 2 ),fi(kx) is satisfied. It corresponds to a 
Diophantine equation but it does not possess the nice properties we used to solve Lemma 8. 
We know moreover that (^1,^2) should be even. A direct search produces the picture given 
in Figure 2. 

Theorem 13. The third order integrability condition for E± is never satisfied except for 




E N 



(23) 



{kt,k 2 ) = (2,2). 
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FIG. 2. Each dot corresponds to a possible even (ki,k 2 ). The solutions seems to be unbounded, 
and the set is probably Zarisky dense. Thus, no practical algebraic information can be extracted 
from this constraint. There are infinitely many solutions (because the diagonal part reduces to a 
Pell equation) with simultaneously arbitrarily high k\,k 2 . So here Theorem 2 is useless without 
Theorem 3. 

Proof. Recall that the parameters (hi, k 2 ) need to be both even for a potential £4 to be inte- 
grable at order 2 near all Darboux points. We begin by solving Qk 2 M {fifa), hi^s), fzik 2 )) = 
in k\. This is a polynomial of degree 4 in k± and as a polynomial, its Galois group is -D4. 
This allows us to write the solution in a relatively simple form 



where Fi,F 2 G Q(/i, /2, fa, k 2 ). Moreover, ki, k 2 are even integers. Let us prove that in fact, 
for even k 2 > 200, the expression 



is always complex for all possible valuations of the square roots. To have real values, we 
need that Fi(k 2 ) + wF 2 (k 2 ) be positive for at least one valuation of the square root. Let us 
begin by proving that w never vanishes. The function w 2 is a polynomial in Q[/i, f 2 , ^3, k 2 ]. 




(24) 



9(k 2 + 2f{k 2 - l) 2 h{k 2 )f 2 {k 2 ) - 6(k 2 + 3)(fc 2 - 2)f 2 {k 2 )h{k 2 ) + 36f 3 (k 2 ) 2 
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Thanks to Theorem 3, we can express /i, /2, /3 in k 2 with controlled relative error. We check 
that the amplication of the error is small after summation of all terms (here it is less than 
1 + 10~ 3 ) and that the approximated expressions never vanish. Now we need to prove that 

Fi(fc 2 ) + wF 2 (k 2 ) < and F x {k 2 ) - wF 2 (k 2 ) < 0. 

We first prove that F 2 {k 2 ) and Fi{k 2 ) (which are in Q[/i, f 2 , fa, k 2 ] of degree 3,4 in fa 
respectively) are always negative. Then we just have to prove that 

4^ > 1 <=> 4^ > 1 <=> w - > <>■ 

wF 2 (k 2 ) w 1 F 2 [k 2 y 
The last expression is in Q[/i, fa, fa, k 2 ] (of degree 8 in fa), so we can prove this statement. 
Again we compute the error amplification of the sum and it stays below 1 + 10~ 3 , and the 
error is then still less than 10~ 4 . Eventually, we prove that this approximated expression 
never vanishes and is always positive. For the remaining cases, we use exhaustive testing 
and we find only one solution [k\, k 2 ) = (2, 2). □ 

The case (ki,k 2 ) = (2,2) corresponds to the second case of Theorem 4. It is really 
integrable with a quadratic in momenta additional first integral which is given in 7 page 107 
case (8). 

IX. REMAINING CASES AND CONCLUSION 

The remaining cases are the ones which do not posess a non-degenerate Darboux point. 

Theorem 14. Consider the set of potentials V given by (6) and suppose that V does not 
possess a non- degenerate Darboux point c. IfV is meromorphically integrable, then V belongs 
to one of the families 

V= l -(a + be ie ), V= l -( a + he ^ i 

V= l -(a + be™), V = ^ a + be ^\ 

with aeC,beC*. 

Proof. First let us suppose that V does not possess any Darboux point c. This means that 
the function 

U{6) = a + be id + ce 2W + de m 
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does not possess any critical point. The only possibility is that U{9) = F(e 10 ) with F(z) = 
a + bz 11 , b 7^ 0. This corresponds to the three first cases of Theorem 14. Now suppose 
there exists one Darboux point c but degenerate. After rotation, we can suppose that the 
Darboux point corresponds to 9 = 0. We have moreover the integrability constraint that 
U"(0) = 0. This gives the potential 

After rotation, this corresponds to the fourth case of Theorem 14. □ 

Remark that in contrary to previous sections, we use the hypothesis that the first integrals 
are meromorphic on C 2 x 5, including r = 0. If we consider first integrals only meromorphic 
on C 2 x5*, the integrability constraint U"(0) = does not hold anymore and two additional 
cases appear in Theorem 4 

v = I ( a + heW ) 2 > v = I { a + bei6 f { a ~ 2heW ) 

The family V = £ (a + be l6 ^j is integrable as given in 7 , but with an additional first integral 
only meromorphic on C 2 x S*. For the other ones, the integrability status is still unknown. 
Let us remark now on the open cases. After rotation and dilatation, these cases correspond 
in fact to a finite number of potentials which are the following (the last two cases being open 
only for meromorphic first integrals on C 2 x S*): 

y = r -i e 2^ y = T -if e 2ie -i) v = r~ 1 ( e m - l) , 

V = r" 1 (e id - if , V = r- 1 (e ie - if , V = r" 1 (e id - if (2e ie + l) . 

We cannot study these potentials because we do not have a particular solution to study, 
or a sufficiently non-degenerate one (studying degenerate Darboux points with higher vari- 
ational method is in fact useless and does not give any additional integrability condition). 
This is of course the main weakness of the Morales-Ramis theory. This is not due to the 
difficulty of applying the Morales-Ramis theory as we treat it in this article, but much more 
a fundamental limitation that seems hard to overcome. One approach could consist in look- 
ing for special algebraic orbits of these systems using a direct search (following Hietarinta 7 ). 
This is not successful for all these potentials. 

To conclude, let us remark that our holonomic approach to higher variational methods 
is very general, and in no way limited to this example. This could work at least for all 
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problems about integrability of homogeneous potentials, as it allows to compute various 
higher integrability conditions of any fixed order. This is linked to the fact that the first order 
variational equation of a natural Hamiltonian system often corresponds to a spectral problem 
of a second order differential operator, which generates P-finite sequences of functions, which 
in turn appear in the study of higher variational equations. We could also wonder if these 
arbitrary high eigenvalues are really possible, and if this work is only conceptual and in 
practice useless. Indeed, very high eigenvalues should correspond to very high degree first 
integrals, and counting the number of conditions and number of free parameters for the 
existence of such high degree first integrals strongly suggests they do not exist. But this 
intuition is wrong, as Andrzej J. Maciejewski, Maria Przybylska found quite recently such 
an example in dimension 3. This is probably linked to the fact that most of integrable cases 
come from ultra-degenerate cases, as in our analysis: the generic case E4 contains only one 
possibility, and when we look at the third order integrability condition, it seems really to 
be a miracle that this condition could ever be satisfied. On the contrary, the cases without 
Darboux points contain lots of integrable potentials. 
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